Generalized BBV Models for Weighted Complex Networks 
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We will introduce two evolving models that characterize weighted complex networks. Though 
the microscopic dynamics are different, these models are found to bear a similar mathematical 
framework, and hence exhibit some common behaviors, for example, the power-law distributions 
and evolution of degree, weight and strength. We also study the nontrivial clustering coefficient 
C and tunable degree assortativity coefficient r, depending on specific parameters. Most results 
are supported by present empirical evidences, and may provide us with a better description of the 
hierarchies and organizational architecture of weighted networks. Our models have been inspired by 
the weighted network model proposed by Alain Barrat et al. (BBV for short), and can be considered 
as a meaningful development of their original work. 

PACS numbers: 



I. INTRODUCTION 



The recent few years have witnessed a great develop- 
ment in physics community to explore and characterize 
the underlying laws of complex networks, including as 
diverse as the Internet the World-Wide Web Q|, the 
scientific collaboration networks (SCN) 0, , and world- 
wide airport networks (WAN) ■ Many empirical mea- 
surements have uncovered some general scale-free prop- 
erties of those real systems, which motivated a wealth 
of theoretical efforts devoted to the characterization and 
modelling of them. Since Barabasi and Albert introduced 
their seminal BA model, most efforts have been con- 
tributed to study the network topological properties • 
However, networks as well known are far from boolean 
structures, and the purely topological representation of 
them will miss important attributes often encountered 
in real world. For instance, the amount of traffic char- 
acterizing the connections of communication systems or 
large transport infrastructure is fundamental for a full 
description of them || . More recent years, the availabil- 
ity of more complete empirical data and higher computa- 
tion ability permit scientists to consider the variation of 
the connection strengths that indeed contain the physical 
features of many real graphs. Weighted networks can be 
described by a weighted adjacency N x N matrix, with 
element Wij denoting the weight on the edge connecting 
vertices i and j. As a note, this paper will only consider 
undirected graphs where weights are symmetric. As con- 
firmed by measurements, complex networks often exhibit 
a scale-free degree distribution P{k) ~ fc~ 7 with 2< 7 <3 
[EEl- Interestingly, the weight distribution P(w) is also 
found to be heavy tailed in some real systems As 
a generalization of connectivity ki, the vertex strength 



is defined as s, 
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j, where T(i) denotes the 



set of i's neighbors. This quantity is a natural measure 
of the importance or centrality of a vertex in the net- 
work. For instance, the strength in WAN provides the 
actual traffic going through a vertex and is obvious mea- 
sure of the size and importance of each airport. For the 
SCN, the strength is a measure of scientific productivity 
since it is equal to the total number of publications of any 
given scientist. Empirical evidence indicates that in most 
cases the strength distribution has a fat tail [(|, similar 
to that of degree distribution. Highly correlated with the 
degree, the strength usually displays scale-free property 
s ~ k 13 with P > 1 0,^3- Driven by new empirical find- 
ings, Alain Barrat et al. have presented a simple model 
(BBV for short) that integrates the topology and weight 
dynamical evolution to study the dynamical evolution 
of weighted networks [l^. An obvious virtue of their 
model is its general simplicity in both mechanisms and 
mathematics. Thus it can be used as a starting point 
for further generalizations. It successfully yields scale- 
free properties of the degree, weight and strength, just 
depending on one parameter 5bbv that controls the lo- 
cal dynamics between topology and weights. Inspired by 
BBV's work, a class of evolving models will be presented 
in this paper to describe and study specific weighted net- 
works. This paper is organized as follows: In Section 
II, we will introduce a traffic-driven model to mimic the 
weighted technological networks. Analytical calculations 
are in consistent with numerical results. In Section III, 
a neighbor-connected model is proposed to study social 
networks of collaboration, with the comparison of simu- 
lations and theoretical prediction as well. At the end of 
each section, we discuss the differences between the BBV 
model and ours, from the microscopic mechanisms to ob- 
served macroscopic properties. We conclude our paper 
by a brief review and outlook in Section IV. 
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II. MODEL A 



The Traffic-Driven Model for Technological 
Networks 



The network provides the substrate on which numerous 
dynamical processes occur. Technology networks pro- 
vide large empirical database that simultaneously cap- 
tures the topology and the traffic dynamics taking place 
on it. We argue that traffic and its dynamics is a key role 
for the understanding of technological networks. For In- 
ternet, the information flow between routers (nodes) can 
be represented by the corresponding edge weight. The to- 
tal (incoming and outgoing) information that each router 
deals with can be denoted by the node strength, which 
also represents the importance or load of given router. 
Our traffic-driven model starts from an initial configura- 
tion of Nq vertices fully connected by links with assigned 
weight wq = 1. The model is defined on two coupled 
mechanisms: the topological growth and the increasing 
traffic dynamics: 

(i) Topological Growth. At each time step, a new ver- 
tex is added with m edges connected to m previously 
existing vertices (we hence require No > m), choosing 
preferentially nodes with large strength: i.e. a node i is 
chosen by the new according to the strength preferential 
probability: 



n 



Si 



new — 



(i) 



The weight of each new edge is also fixed to Wo = 1 • This 
strength preferential mechanism have simple physical and 
realistic interpretations in Ref. 0, 0] . 

(ii) Traffic Dynamics. From the start of the net- 
work growing, traffic in all the sites are supposed to 
constantly increase, with probability proportional to the 
node strength Sij ^ fc Sk per step. We assume the grow- 
ing speed of the network's total traffic as a discrete con- 
stant W (each unit can be considered as an information 
packet in the case of Internet). Then in statistic sense 
the newly created traffic in site i per step is 



AW, = W 



Efe s fe 



(2) 



These newly-added packets will be sent out from i to their 
separate destinations. Our model does not care their spe- 
cific destinations or delivering paths, but simply suppose 
that each new packet preferentially takes the route-way 
with larger edge weight (data bandwidth of links), i.e. 
with the probability Wij/si, and it hence will increase the 
traffic (strength) in the corresponding neighbor j of node 
i. It is a plausible mechanism in many real- world webs. 
For instance, in the case of the airport networks, the 
potential passenger traffic in larger airports (with larger 
strength, often located in important cities) will be usually 
greater than that in smaller airports, and busy airlines 
often get busier in development. For Internet, routers 



that have larger traffic handling capabilities are respon- 
sible to deal with more information packets. Also, the 
route-ways with broader data bandwidth will get busier. 
Admittedly, this "busy get busier" scenario is intuitive 
in physics, though perhaps not strict in mathematics. 



B. Analytical Results vs. Numerical Simulations 

The model time is measured with respect to the num- 
ber of nodes added to the graph, i.e. t = N — Nq, and the 
natural time scale of the model dynamics is the network 
size N. In response to the demand of increasing traffic, 
the systems must expand in topology. With a given size, 
one technological network assumably has a certain abil- 
ity to handle certain traffic load. Therefore, it could be 
reasonable to suppose for simplicity that the total weight 
on the networks increases synchronously by the natural 
time scale. That is why we assume W as a constant. 
This assumption also bring us the convenience of analyt- 
ical discussion p|. By using the continuous approxima- 
tion, we can treat fc, w, s and the time t as continuous 
variables. The time evolution of the weights Wij can be 
computed analytically as follows: 



dwij 
~dT 
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(3) 



The term AWiWij / Si represents the contribution to 
weight w^ from site i. Considering 



J^s fc w (2W + 2m)t, 



one can rewrite the above evolution equation as: 



dwij W 



alt W + m t 



(4) 



(5) 



The link (i,j) is created at tij = max(i,j) with initial 
condition to, j (f y ) = 1, so that 



(i) 



, x W/(W+m) 



(6) 



Further, we can obtain the evolution equations for Si(t) 
and ki(t): 
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These equations can be readily integrated with initial 
conditions fc,(£j) = Si(ti) = m, yielding 



Si(t) 



t 

i 

2W 



2W + m 



(9) 
(10) 



The strength and degree of vertices are thus related by 
the following expression 



2W 



2W. 



(11) 



In order to check the analytical predictions, we performed 
numerical simulations of networks created by the present 
model with various values of W and minimum degree m. 
In Fig. 1, we report the average strength Si of vertices 
with connectivity ki and confirm the validity of Eq. (11). 

The time ti — i when the node i enters the system is 
uniformly distributed in [0, t] and the strength probabil- 
ity distribution can be written as 



P(s,t) = 



1 



t + N 



5(s - s^t^dti 



(12) 



where S(x) is the Dirac delta function. Using equation 
Si(t) ~ (t/i) a obtained from Eq. (9), one obtains in 
the infinite size limit t — > oo the power-law distribution 
P(s) ~ s^a (as shown in Fig. 2) with 



T A = 1 + - = 2 
a 



2W 



(13) 



Obviously, when W = the model is topologically equiv- 
alent to the BA network and the value j A — 3 is recov- 
ered. For larger values of W, the distribution is gradually 
broader with "f A — > 2 when W — > oo. Since s and k are 
proportional, one can expect the same behavior of degree 
distribution P(k) ~ . Analogously, the weight dis- 
tribution can be calculated yielding the scale-free prop- 
erty P(w) ~ w~ a A 1 with the exponent a A = 2 + m/W 
(See Fig. 3). 



C. Clustering and Correlations 




FIG. 1: The average strength Si of the nodes with connec- 
tivity ki for different W (log- log scale). Linear data fittings 
all give slope 1.00 ± 0.02, demonstrating the predicted linear 
correlation between strength and degree. 




FIG. 2: Cumulative strength probability distribution P c (s) 
with various W. Data agree well with the power-law form 
s iA . The inset reports j A from data fitting (filled circles), in 
comparison with the theoretical prediction 7 A = 2 + 2 w+m = 
5 r(line), averaged over 20 independent networks of size 



2+: 



2W+5 1 

N=7000. 



vertices' triples in the network. Further information can 
be gathered by inspecting the average clustering coeffi- 
cient C(j) restricted to classes of vertices with degree k: 



A complete characterization of the network structure 
must take into account the level of clustering and degree 
correlations present in the network. Information on the 
local connectedness is provided by the clustering coeffi- 
cient Ci defined for any vertex i as the fraction of con- 
nected neighbors of i. The average clustering coefficient 
C = TV" 1 J2i c» thus expresses the statistical level of co- 
hesiveness measuring the global density of interconnected 



v ; i/ki—k 

In many networks, C(k) exhibits a power-law decay as 
a function of k, indicating that low-degree nodes gen- 
erally belong to well interconnected communities (high 
clustering coefficient) while high-degree sites are linked 
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FIG. 3: Weight probability distribution P(w) with various 
values of W. Data are consistent with a power-law behavior 
w~ aA . In the inset we give the value of a A obtained by data 
fitting (filled circles), together with the analytical expression 
a A — 2 + m/W — 2 + 5/W (line). The data are averaged over 
20 independent realizations of network size N=7000. 



FIG. 4: The scaling of C(k) with k for various W. The 
data are averaged over 20 independent runs of network size 
N=7000. The inset shows the average clustering coefficient 
C, depending on increasing W. 



to many nodes that may belong to different groups which 
are not directly connected (small clustering coefficient). 
This is generally the signature of a nontrivial architecture 
in which hubs (high degree vertices) play a distinct role 
in the network. Correlations can be probed by inspecting 
the average degree of the nearest neighbors of a vertex i, 
that is, k nn ,i = k^ 1 Y^j kj. Averaging this quantity over 
nodes with the same degree k leads to a convenient mea- 
sure to investigate the behavior of the degree correlation 
function 



k n n (fc) 



1 



NP{k) 



J2 k nn ^ = J2k'P(k'\k), (15) 



i/ki—k 



If degrees of neighboring vertices are uncorrelated, 
P(k'\k) is only a function of k' and thus k nn (k) is a con- 
stant. When correlations are present, two main classes 
of possible correlations have been identified: assortative 
behavior if k nn (k) increases with k, which indicates that 
large degree vertices are preferentially connected with 
other large degree vertices, and dis assortative if k nn (k) 
decreases with k. The above quantities provide clear sig- 
nals of a structural organization of networks in which 
different degree classes show different properties in the 
local connectivity structure. Almost all the social net- 
works empirically studied show assortative mixing pat- 
tern, while all others including technological and biolog- 
ical networks are disassortative. The origin of this dif- 
ference is not understood yet. To describe these types 
of mixing, Newman further proposed some simpler mea- 
sures, which is called assortativity coefficients [l4|. In 
this paper, we will also use the formula proposed by New- 




FIG. 5: Average connectivity k nn (k) of the nearest neighbors 
of a node depending on its connectivity k for different W . The 
data are averaged over 20 network realizations of N=7000. 
The assortativity coefficient r depending on W is shown in 
the inset. 



man UJ, 



(16) 



where ji , k{ are the degrees of vertices at the ends of the 
ith edges, with i = 1, ...,M (M is the total number of 
edges in the observed graph) . Simply, r > means assor- 
tative mixing, while r < implies disassortative pattern. 

In order to inspect the above properties we perform 
simulations of graphs generated by the model for dif- 
ferent values of W, fixing TV = 7000 and m = 5. In 
the case of clustering, the model also exhibits properties 
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which are depending on the basic parameter W. More 
precisely, for small W, the clustering coefficient of the 
network is small and C(k) is flat. As W increases how- 
ever, the global clustering coefficient increases and C(k) 
becomes a power-law decay similar to real network data 
15]. Fig. 4 shows that the increase in clustering is de- 
termined by low-degree vertices. The average clustering 
coefficient C is found to be a function of W, as shown in 
its inset. This numerical result obviously demonstrates 
the important effect of traffic on the hierarchical struc- 
ture of technological networks. Analogous properties are 
obtained for the degree correlation spectrum. For small 
W, the average nearest neighbor degree k nn {k) is quite 
flat as in the BA model. The disassortative character 
emerges as W increases and gives rise to a power law 
behavior of k nn (k) ~ k~ n (Fig. 5). The assortativity 
coefficient r versus W, as reported in the inset of Fig. 
5, demonstrates the tunable disassortative property of 
this model, which is supported by measurements in real 
technological networks The qualitative explanations 
of the correlations and clustering spectrum can be found 
in [J^, and their theoretical analysis appears in [l7| . In 
sum, all the simulation results for clustering and degree 
correlation, as empirically observed, imply us that the 
increasing traffic may be the driven force to shape the 
hierarchical and organizational structure of real techno- 
logical networks. 



D. Comparison with the BBV model 

One may notice that if the parameter W is replaced 
by W = m5 A , then the mathematical framework of our 
traffic-driven model is equivalent to that of BBV's [lfij . 
though the specific weights' dynamics are quite different 
between the two. By comparison with S BBV (that is, the 
fraction of weight which is locally "induced" by the new 
edge onto the neighboring others), the parameter 5 A in 
our model, with macroscopic perspectives, is the ratio 
of the total weight increment on the expanding struc- 
ture (W) to the number of newly-established links (m) 
at each time step. It is an important measure for the 
relative growing speed of traffic vs. topology, and con- 
trols a series of the network scale-free properties. Thus, 
there is an obvious difference between BBV model and 
ours: the former is based on the local rearrangement of 
weights induced by newly added links, while the latter 
is built upon the global traffic growth and the redistri- 
bution of weights according to the local nature of net- 
work. Noticeably, the weight dynamical evolution of the 
BBV model is triggered only by newly added vertices, 
hardly resulting in satisfying interpretations to collabo- 
ration networks or the airport systems. For these two 
cases, even if the size of networks keeps invariant, co- 
authored papers will still come out and airports can be- 
come more crowded as well. In contrast, this difficulty for 
practical explanations naturally disappear in our model, 
due to its global weight dynamics. Above all, the traffic- 



driven model here, without loss of simplicity and prac- 
tical senses, has generalized BBV's work and narrowed 
its applicable scope to technological networks. Based on 
it, more complicated variations of the microscopic rules 
may be implemented to better mimic technological net- 
works. Especially worth remarking is that the present 
empirical studies on airport networks and Internet in- 
dicate the nonlinear degree-strength correlation s ~ k 13 
with (3 > 1. In Ref. [ly, Barrat et al. proposed the 
heterogeneous coupling mechanism to obtain this prop- 
erty. This is not difficult to introduce into the present 
traffic-driven version. Moreover, we find that the non- 
trivial weight-topology correlation can also emerge from 
the accelerating growth of traffic weight 0] or from the 
accelerating creation and reinforcement of internal edges 
00 



III. MODEL B 

A. The Neighbor-Connected Model for Social 
Networks 

Social networks are a paradigm of the complexity of 
human interactions, which have also attracted a great 
deal of attention within the statistical physics commu- 
nity. The study of social networks has been traditionally 
hindered by the small size of the networks considered and 
the difficulties in the process of data collection (usually 
from questionnaires or interviews). More recently, how- 
ever, the increasing availability of large databases has 
allowed scientists to study a particular class of social 
networks, the so-called collaboration networks. The co- 
authorship network of scientists represents a prototype 
of complex evolving networks, which can be defined in a 
clear way. Their large size has allowed researchers to get 
a reliable statistical description of their topological and 
weight properties, and hence reach reliable conclusions of 
the network structure and dynamics. In weighted social 
networks, the edge weight between a pair of nodes can 
represent the tightness of their connection. The larger 
weight indicates the more frequency of interaction; e.g. 
the number of co-authored papers between two scientists, 
the frequency of telephone or email contacts between two 
acquaintances, etc. In addition, it is more probable that 
two vertices with a common neighbor get connected than 
two vertices chosen at random. Clearly, this property 
leads to a large average clustering coefficient since it in- 
creases the number of connections between the neighbors 
of a vertex. This is already observed in a model proposed 
by Davidsen, Ebel and Bornhodt [lflj . 

Our neighbor-connected model starts from an initial 
graph of No vertices, fully connected by links with as- 
signed weight Wq = 1. Its evolution then is simply based 
on the dynamics of connecting nearest-neighbors: At each 
time step, a new vertex n is added with one primary link 
and m secondary links, which connect with m+1 existing 
vertices (the initial network configuration hence requires 
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No > m + 1). Actually, the newly-built connections are 
not independent, but related with each other. The pri- 
mary link (wq = 1) first preferentially attaches to an old 
node with large strength; i.e. a node i is connected by 
the primary link with probability 



(17) 



Then, the to secondary connections (assigned wq = 1 
each) are preferentially built between the new vertex and 
to neighbors of node i, with the weight preferential prob- 
ability 



(18) 




where j € T(z). After building the primary link (n, i), 
the creation of every secondary link (n,j) is assumed to 
introduce variations of network weights. For the sake 
of simplicity, we limit ourselves to the case where the 
introduction of a primary link on node i will trigger only 
local rearrangement of weights on the existing neighbors 
j E r(i), according to the rule 



FIG. 6: Weight probability distribution with m=3. Data are 
consistent with a power-law behavior P(w) ~ w~ a B . In the 
inset we give the value of a B obtained by data fitting (filled 
circles), together with the analytical form Eq. (22) (line). 
The data are averaged over 100 independent realizations of 
network size N=7000. 



W; 



(19) 



In general, S depends on the local dynamics and can be 
a function of different parameters such as the weight Wij , 
the degree or the strength of i, etc. In the following, 
we will simply focus on the case that <5 = const. After 
the weights have been updated, the evolving process is 
iterated by introducing a new vertex until the desired size 
of the network is reached. 

The above mechanisms have simple physical and real- 
istic interpretations. Once a fresh member joins a social 
community, he will be introduced to his neighbors or take 
initiatives to interact with them. Then, his social connec- 
tions will soon be built within his neighborhood. It is rea- 
sonable that the entering of this new member can trigger 
the local variation of connections. Take the SCN for ex- 
ample, a scientist joining a research group will often col- 
laborate with both the group director (primary link) and 
the other group members (secondary links). The scien- 
tist's affiliation can naturally boost the research produc- 
tivity of the group and also, enhance the collaborations of 
other members. The above scenario may be a best inter- 
pretation of the origin of our model parameter S. We use 
S to control the effect of the newly-introduced member 
on the weights of connections among the local neighbors. 
In the proceeding section, we will see the model also gives 
a wealth of scale-free behaviors depending on this basic 
parameter. 



B. Analytical and Numerical Results 

Along the analytical lines used in Section II, one can 
also calculate the time evolution of the strength, weight 



and degree, and hence calculated their scale-free distri- 
butions as follows: 



dwij msi Wij msj toy w i3 

= — -5 + — -o — 2md- 

dt 
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(20) 



By noticing J2k s k ~ 2(m + l)t + 2m5t = 2(1 + m+mS)t, 
we have 



dwi 3 
~dt 



m8 



1 + m + mS t 



(21) 



and so that Wij(t) = (t/tij) , which indicates the power- 
law distribution of weight P(w) ~ w~~ a B with exponent 



a- = 1 



= 2 



1 + m 



(22) 



Further, the evolution equations for Sj and ki are ob- 
tained 



dsj 
dt 



Edwa 



dt 
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FIG. 7: Strength probability distribution with m=3. Data 
agree well with the power-law form P(s) ~ s 'b . The inset 
reports 7 B from data fitting (filled circles) averaged over 100 
independent networks of size N=7000, in comparison with the 
theoretical prediction Eq. (26) (line). 



FIG. 8: The average strength Sj of vertices with connectivity 
ki for different 5 with m=3 (log- log scale). Linear data fittings 
all give slope 1.00 ± 0.02, confirming the prediction /3 = 1. 



Integrating with the initial conditions ki(t = i) 
i) = m + 1, we have 



i(t) = (m + 1) 



ki(t) = (m + 1) 



t 



Si(t) + 2m5 
1 + to + 2m6 ' 



(25) 



which determine the scale-free distributions of both 
strength and degree, with the same power-law exponent 



7 B =l + i = 2 



1 



1 + to + 2mS 



(26) 



We also performed numerical simulations of networks 
generated by the model with various values of S and min- 
imum degree to. As one can see, Fig. 6 and 7 recover 
the theoretical predictions for scale-free distributions of 
weight and strength, and Fig. 8 validates the linear 
strength-degree correlation. It is worth stressing that 
the empirical evidence in co-authorship networks just in- 
dicates the linear correlation s ~ k la. 



C. Discussions 

There are two important differences between the BBV 
model and our neighbor-connected one, though the lat- 
ter can be regarded as an interesting variation of the 
former. First, in the BBV model the evolution and dis- 
tributions of such quantities as strength, weight and de- 
gree are simply depending on the parameter S BBV which 
controls the coupling between local topology and weights. 
In our model, however, the evolution and distributions of 
those quantities are controlled by two parameters (S and 




FIG. 9: The scaling of C(k) with k for various S with m=3. 
Their fitted power-law exponents all give 0.94 (dashed line). 
The data are averaged over 100 independent runs of network 
size N=7000. The inset shows the average clustering coeffi- 
cient C depending on increasing 5. 



to, as reflected in theoretical part), which together deter- 
mine the effect of the new node on the local weights and 
topology. Though we have fixed m in most simulations, 
the role of parameter to should not be ignored. As Fig. 9 
reports, the curvature of C(k) is not sensitive to the vari- 
ations of 6, but it nontrivially depends on to as shown in 
Fig. 10. Second, given the same minimum degree, the av- 
erage clustering coefficient C of our neighbor-connected 
model (see the inset of Fig. 9) can be much larger than 
that of BBV's model or of our Model A, because the sec- 
ondary links in Model B considerably increase the density 
of triangles within the system. Compared with the orig- 
inal model, this larger clustering demonstrate an impor- 
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FIG. 10: The scaling of C(k) with k for various m with 
5 = 1. The data are averaged over 100 independent runs of 
network size N=7000. The inset shows the average clustering 
coefficient C depending on increasing m. 



the model means larger minimum degree. When a new 
node arrives, the number of triangles in the system will 
increase by m. But for a smallest-degree node i, its clus- 
tering is Cj(m) = 2N A (i)/m(m + 1) where N A (i) denotes 
the number of connected neighbors of node i. Therefore, 
increasing m, we will decrease Cj more greatly, resulting 
in the decaying of C. Admittedly, there still exists a com- 
mon point which leads to the restriction of Model B to 
mimic and interpret social networks. The degree corre- 
lations in both models are negative (see the inset of Fig. 
11), indicating their disassortative mixing patten that is 
opposite in social networks. Is it possible to find out a 
unified minimum model to characterize both the assorta- 
tive and disassortative networks? This question is very 
challenging and appears among the leading ones in front 
of network researchers Q. Our recent studies may shed 
some new light on this tackling problem 20]. Anyway, 
the present neighbor-connected model (though disassor- 
tative) has maintained the simplicity of the BBV model 
and appears as a more specific version for weighted so- 
cial networks, considering its larger clustering and clearer 
hierarchical structure. 




FIG. 11: Average connectivity k nn (k) of the nearest neigh- 
bors of a node depending on its connectivity k for different 
8 with m=3. The data are averaged over 100 network real- 
izations of N=7000. The assortativity coefficient r versus 5 is 
also reported in the inset. 

tant advantage of Model B in modelling the small-world 
property of real complex networks. One question arises 
from the inset of Fig. 10: at first sight, it may be sur- 
prising to see that C greatly decreases when increasing 
the secondary linking number m. Actually, larger m in 



IV. REVIEW AND OUTLOOK 

In this paper, we have presented and studied two evolv- 
ing models for weighted complex networks. These two 
models intend to mimic technological networks and so- 
cial graphs respectively, and can be regarded as a mean- 
ingful development of the original BBV model. Though 
their specific evolution dynamics are different, all of them 
are found to bear similar mathematical structures and 
hence exhibit several common behaviors, e.g. the power- 
law distributions and evolution of degree, weight and 
strength. In each case, we also studied the nontrivial 
clustering coefficient C and tunable degree assortativ- 
ity r as well as their degree-dependent correlations. In 
such context, we compared our generalized models with 
the original one, and got several interesting conclusions, 
which may provide us with a better understanding of the 
hierarchical and organizational architecture of weighted 
networks. For all the above reasons, we would like to clas- 
sify these models into a class, within which the BBV's 
is the first one and may be the simplest one. It must 
be admitted that this class of models do not take into 
account the internal connections during the network evo- 
lution, which is yet beyond the scope of this paper. Still, 
we hope that our generalized work can make this model 
family more diversified and help bring it some new sights. 
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